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Abstract 



We compute the Lie symmetry algebra of the generalized Davey-Stewartson 
(GDS) equations and show that under certain conditions imposed on param- 
eters in the system it is infinite-dimensional and isomorphic to that of the 
standard integrable Davey-Stewartson equations which is known to have a very 
specific Kac-Moody-Virasoro loop algebra structure. We discuss how the Vi- 
rasoro part of this symmetry algebra can be used to construct new solutions, 
which are of vital importance in demonstrating existence of blow-up profiles, 
from known ones using Lie subgroup of transformations generated by three- 
dimensional subalgebras, namely sl(2, R). We further discuss integrability as- 
pects of GDS equations. 



1 Introduction 

A system of nonlinear partial differential equations in 2+1 dimensions as a model 
of wave propagation in a bulk medium composed of an elastic material with couple 
stresses has recently been derived in Ref. PQ, namely 
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with the condition (A — l)(mi — 1112) = n 2 . Here if)(t,x,y) is a complex function, 
w(t,x,y) and <p(t,x,y) are real functions and S, n, mi, m-z, A, x, 7 are real constants. 
The authors of PP showed that if the parameters are related by 



then (jl.lj) can be reduced to the standard Davey-Stewartson (DS) equations (in gen- 
eral not integrable) by a non-invertible point transformation of dependent variables. 
Therefore, they called the generalized Davey-Stewartson (GDS) equations. Be- 
low, we justify this naming from a group-theoretical point of view. Also, in jl, 
some travelling type solutions of (jl.lj) in terms of elementary and elliptic functions 
are obtained. Based on some physically obvious Noetherian symmetries (time-space 
translations and constant change of phase), global existence and nonexistence results 
are given in 0. In another recent work under some constraints on the physi- 
cal parameters, the so-called hyperbolic-elliptic-elliptic case of the system (jl.lj) (In 
[3] the system is classified into different types according to the signs of parameters 
(S, mi, m 2 , A)) was shown to admit singular solutions that blow up in a finite time. 
To do this, inspired by the (pseudo) conformal invariance of DS system, they used 
the fact that time-dependent SL(2,M) invariant solutions can be generated from sta- 
tionary radial solutions for an appropriate choice of coefficients. 

The purpose of this article is to study GDS equations from group theoretical point 
of view. For our purposes, we find it more convenient to consider the differentiated 
form of (jl.lj) . Thus, differentiating the last two equations of (jl.lj) with respect to x 
and y, respectively and then making the substitution w x — > w, (j> y — > 4> and rewriting 
the corresponding system in a real form by separating ip = u + iv into real and 
imaginary parts, we obtain a system of four real partial differential equations 



In the sequel, we shall call (jl.3j) the GDS equations. 

The main result of the paper is to show that, when some conditions on physical 
parameters 5, n, mi, m 2 , A are imposed, the Lie algebra of the symmetry group of the 
GDS system has a Kac- Moody- Virasoro (KMV) loop structure which is shared by the 
symmetry algebras of all known integrable equations in 2+1 dimensions such as the 
Kadomtsev-Petviashvilli (KP) equation |U |3] and the usual integrable DS equations 
jHj. The corresponding special case is candidate for being integrable. Moreover, we 



n = 1 — A = mi — m 2 



(1.2) 



u t + Sv xx + v yy = xv{u 2 + v 2 ) + ^v(w + (j)), 
-v t + 5u xx + u yy = x u (u 2 + v 2 ) + ju(w + <p), 
w xx + n<j) xx + m 2 Wyy = 2(u 2 x + uu xx + v 2 x + vv xx ), 

UWyy + \(f) XX + m X (j)yy = 2^ + UUyy + vl + VVyy) . 



(1.3) 
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show that this algebra is isomorphic to that of DS equations [Bj 

ilpt + hlpxx + 4>yy = $2\4>\ 2 ij + Wlp 
ExW xx + W yy = S 2 (\i>\ 2 )yy, 

with 5i — ±1, 82 — ±1. In the special case when 8\ + E\ — 0, this system which is 
one of the rare systems in more than 1+1 dimensions for which the Cauchy initial 
value problem is solvable by the inverse spectral transform (1ST) technique becomes 
completely integrable. The Lie algebra of the symmetry group of the integrable DS 
system is referred to as the DS algebra. This isomorphism (a necessary condition 
for two different systems to be transformable into each other) should motivate one 
to look for point transformations taking the Lie algebras into each other. We expect 
such transformations to transform the systems into each other as well. 

In Section 2 we compute the Lie symmetry algebra of (jl.3j) and identify its struc- 
ture. In particular, we show that for special choice of parameters it is a centerless 
Kac-Moody-Virasoro algebra. 



2 The symmetry group of the GDS equations and 
structure of its Lie algebra 

We apply the standard infinitesimal procedure [7j to find the symmetry algebra L 
and hence the symmetry group G of (jl.3j) . We write the GDS equations as a system 
Aj(£, x, y, u, v, w, <fi) = 0, i = 1, 2, 3, 4. A general element of the algebra is represented 
by a vector field 



V = rd t + (d x + r]d y + (fA + tp 2 d v + ip 3 d u 



(2.1) 



where the coefficients T,£,,r),(fi,i = 1,2,3,4 are functions of t,x,y,u,v,w,(f). Ac- 
cording to the general theory for symmetries of differential equations, to find these 
functions we prolong the vector field (j2.1j) to second order derivatives and require 
that the second prolonged vector field annihilates A 4 on the solution manifold of the 
system, namely 



pr 



{2) V(Ai(t,x,y,u,v,w,<p)) 



0. 



A,;=0 



1,2,3,4, 



(2.2) 



where pr^ V is the second prolongation of the vector field V. This condition provides 
us with a quite complicated system of determining equations (a system of linear partial 
differential equations) for the coefficients. This step is entirely algorithmic and is 
implemented on several computer algebra packages like REDUCE, MATHEMATICA, 
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MAPLE (See jS] for a survey of symbolic softwares for symmetry). The final step of 
integrating the determining equations is less algorithmic. Solving these huge number 
of determining equations we find that the general element can be written as 



V = T(f) + X(g) + Y{h) + W(m), (2.3) 



where 



T(f) = f{t)dt + ~f'(t)(xd x + yd y - ud u - vd v - 2wd w - 

X(g) = g(t)d x - ^- § [g'(t)(vd u - ud v ) + 9 -^-{d w + fy)], (2.4) 

= /i(t)^ - - + ^(a w + a^)], 

z 27 
m'(t) 

W(m) = m(t)(vd u - ud v ) + —^(d w + d+). 

2 7 

The functions g{t), h(t), and m(t) are arbitrary functions of class C°°(I), ICR. The 
function f(t) is arbitrary if 

m 2 5 + n + 1 = 0, rrixd + n5 + X = 0, (2.5) 

otherwise f(t) = c^i 2 + Cit + Cq. We mention that these conditions come from the fact 
that two of the determining equations are 

(m 2 S + n + l)f"'(t) = 0, {m x 5 + nd + X)f"(t) = 0, 

whereas the remaining ones are solved without any constraints on g, h and m. 

We mainly focus on the case when f(t) is allowed to be arbitrary. The symmetry 
algebra realized by the vector fields ()2.3j) and ()2.4j) is then infinite-dimensional and 
more important has the structure of a Kac-Moody-Virasoro algebra as we shall see 
below. More interestingly, it is generic among the symmetry algebras of a few 2+1- 
dimensional integrable partial differential equations (the KP equation, the modified 
KP equation, the potential KP equation, the integrable three-wave resonant equations 
and the integrable DS equations). Henceforth, we shall call this the GDS symmetry 
algebra and the corresponding system the GDS system. 

Note that sometimes it is more convenient to use the polar decomposition u + iv = 
Re la so that in (j2.4j) we can write 

ud u + vd v = Rd R , -(vd u - ud v ) = d a . 
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The commutation relations for the GDS algebra are easily obtained as follows: 



[T{fi),T(J 2 )]=T(f 1 &-f 1 f 2 ) 
[T(f),X(g)} = X(fg'- 1 -fg) 

[T{f),Y{h)\=Y{fh'- l -fh) 
[T(f),W(m)] = W(fm') 



(2.6) 



[X( gi ),X(g 2 )] = -j- 6 W(g ig ' 2 - 

[Y(h 1 ),Y(h 2 )] = -^W(h 1 h' 2 - 
[X(g),Y(h)} = [X(g),W(m)] = 



■■ [Y(h), W(m)\ = [W(mi), W(m 2 )} = 0. 



9i92) 



From ()2.6|) we see that the GDS system has a Lie symmetry algebra L isomorphic to 
that of the DS symmetry algebra j^j. Indeed, it allows a Levi decomposition 



is a nilpotent ideal (nilradical). Here, GE- denotes a semi-direct sum. The algebra 
{T(f)} is isomorphic to the Lie algebra corresponding to the Lie group of diffeomor- 
phisms of a real line. 

We remark that a similar isomorphism between the symmetry algebras of a class 
of (integrable) generalized cylindrical KP (GCKP) equation and of the KP equation 
was pointed out in Ref. 

Expanding the arbitrary functions /, g, h and m into Laurent polynomials and 
considering each monomial t n (n not necessarily positive integer) separately, we obtain 
a realization of a KMV algebra without central extension. Here the factor subalgebra 
S is the Virasoro part, the nilpotent subalgebra N is the Kac-Moody part of the 
GDS algebra [TO] - We refer for different realizations of the Virasoro algebras to [11J. 
Furthermore, just as the DS algebra [Hj it can be shown that the GDS algebra with 
(12. 5|) can be imbedded into a Kac-Moody-type loop algebra. 

Theorem 2.1 The system (|1.3|) is invariant under an infinite- dimensional Lie point 
symmetry group, the Lie algebra of which has a Kac-Moody- Virasoro structure iso- 
morphic to the DS algebra if and only if the conditions (12. 5|) hold. 



L = S eiv, 



(2.7) 



where S = {T(f)} is a simple infinite dimensional Lie algebra and 



N={X(g),Y(h),W{m)} 
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Let us mention that the GDS equations are also invariant under a group of discrete 
transformations generated by 



t — > t, x — ► —x, y — > y, -0 — > i/j, w — > w, ~^ 4> 

x->x, y->-y, w->w, -> 

t t, x — > x, y ^ y, "0 — > — 0", w ^ w, (ft —> <fi 

t — > — t, x — > a;, y ^ y, tj; — > -0*, u> — > u>, 0^0. 

The obvious physical symmetries L p of the GDS equations are obtained by re- 
stricting all the functions /, <?, /i and m to be first order polynomials. Indeed, we 
have 

r = r(i) = a = x(i) = a,, p 2 = = a, 
Wo = w(i) = ud u - ud v , 

D = T(t) = td t + \ {xd x + yd y - ud u - vd v - 2wd w - 2<j)d (j> ) 
x 

Bi = X(t) = td x - -^( v du - ud v ), B 2 = Y(t) = td y - -{vd u - ud v ) 
Wt = W(t) = t(vd u - ud v ) + ^-(d w + 8+). 

We see that T,Px,P2 generate translations, D dilations, B\ and B2 Galilei boosts in 
the x and y directions, respectively. Finally, Wo and W\ generate a constant change 
of phase of 0> and a change of phase of ip, linear in t, plus constant shifts in w and 0, 
respectively. 

The generators (|2.9|) form a basis of a eight-dimensional solvable Lie algebra 
L p = {D, T, Pi, P 2 , B\, B 2 , Wo, Wi}. It has a seven-dimensional nilpotent ideal (the 
nilradical) N = {T, P l5 P 2 , B x , B 2 , W , W x }. 

Another finite-dimensional algebra, not contained in L p is obtained by restricting 
f(t) to quadratic polynomials. We obtain T = T(l), D = T(t) as in (J2.9)) . and in 
addition 

C = T(t 2 ) = t 2 d t +tD- (X |/ y \ vd u - ud v ). (2.10) 
The commutation relations are 

[T,D] = T, [T,C] = 2D, [D,C] = C, 
so that we have obtained the algebra sl(2,R) with C generating conformal type of 
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transformations 



t x y 

t - oc ■ y . 

i — pt i — pt i — pt ' 

5/1 ,Nn ~ + , (2.11) 

w; = (1 — pt) 2 w, 4> = (1 — pt) 2 (p, 

where p is the group parameter. Further, composing (|2.11|) with time translations 
generated by T and dilations generated by D we obtain the SL(2,R) group gener- 
ated by actions on the space of independent and dependent variables. It should be 
mentioned that any finite dimensional subalgebra of the Virasoro algebra of 2+1 di- 
mensional integrable equations is isomorphic to sl(2, R) or one of its subalgebras. The 
transformed variables and the new solution in terms of the original ones are given by 
the formulas 

~ c+dt x _ y 7 7 -t 

t— — , x = — , y= — , ad — be = 1 

a + bt a + bt a + bt 

~ , . ._ n r ib(x 2 + Sy 2 ) , , ,~ _ „. 

4}=(a + bt) exp{ 4S , +bt , mt,x,y) ( 2 .12) 

w = (a + bt)~ 2 w(t, x, y) 
(a + bt)- 2 4>(t,x,y). 

Here a,b,c are the group parameters of SL(2,R). These are exactly the formulas 
which played an essential role in the construction of analytic blow-up profiles |3] 
in which the authors made use of stationary radial solutions (ip,w,<p) to generate 
new solutions (time dependent) (if) , w , <ft) of the GDS equations. More generally, the 
elements of the connected part of the full symmetry group of the GDS equations can 
be obtained by integrating the vector fields ([2.3)1 . ([2.41) . We refer the reader to Ref. 
[H] for the general Lie group of transformations of DS algebra. 

Let us now return to the isomorphic GDS and DS symmetry algebras, and trans- 
form the GDS vector fields (j2.ll) by the point transformation q = w + <fi — \ip\ 2 . It is 
easy to see that the component \(d w + d^) transforms to d q , and D — > xd x + yd y — 
ud u — vd v — 2qd q , and the rest remains unaltered, namely the DS symmetry algebra 
is obtained. This means that the functions (if), q) satisfy the DS equations whenever 
(if), w, <p) satisfy the GDS equations, but not vice versa. At this time, it remains open 
whether it is possible to construct an invertible point transformation relating these 
two systems. 

We conclude by making several comments. As is well illustrated by the results of 
this paper, knowing that a nonlinear partial differential equation (or system) admits 
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a KMV algebra as a symmetry algebra can serve as a useful criterion of identifying 
integrable equations. In particular, this fact can be used to pick out an integrable 
equation from a class of generically nonintegrable ones. For instance, for all values 
of parameters not satisfying (J2.5j) . the Virasoro part T(f) of the GDS algebra is 
not present. The first author of the present paper and Winternitz ^2] used the 
same approach to identify all subclasses invariant under a KMV algebra and its 
subalgebras containing up to three arbitrary functions of time from a rather general 
class of KP type equations involving 9 arbitrary functions of one or two variables. On 
the other hand, a classification of all one- and two-dimensional subalgebras of the DS 
algebra into conjugacy classes under the adjoint action of the DS group (including 
the discrete transformations) is performed in [S]. The GDS algebra will have the 
same conjugacy classes of subalgebras as the DS algebra. Depending on which of 
the functions g(t), h(t) and m(t) are nonzero, precisely six conjugacy classes of one- 
dimensional subalgebras exist: 

L ltl = {T(l)}, L? >2 = {X(l)+aY(l)} ; L^h) = {X (1) + Y (h)} , a > 0, 

L M = {y(l)}, L 1>5 = {W(t)}, L lfi = {W(l)}. 

They can be used to reduce the integrable GDS system to integrable one in two 
variables and thus to obtain subgroup invariant solutions. There will be four type of 
reductions since only the first four subgroups corresponding to (£1,1, Lf 2 , L 1(3 (/i), L 14 ) 
will generate actions on the coordinate space (t,x,y). The remaining two (L 15 ,L 16 ) 
generate purely vertical (or gauge) transformations changing phases only and thus 
lead to no reductions. For example, one can show that all the travelling wave solutions 
obtained in [T] can be extracted from those of representative reduced equations by 
applying appropriate symmetry group transformations. We note that these type of 
physically important solutions are invariant under translational subgroups alone. 
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